Strong Dissipation Inhibits Losses and Induces Correlations in Cold Molecular Gases 



OO 

o 
o 

(N 
(N 



5-H 

CD 



I 

O 

o 



> 
o 

o 

00 

o 



X 



N. Syassen/ D. M. Bauer,i M. Lettner,i T. Volz,i'B D. Dietze,i'[l| 
J. J. Garcia- RipoU,^'^ J. I. Cirac/ G. Rempe/ and S. Diirr^ 

^ Max- Planck- InsUtut fur Quantenoptik, Hans-Kopfermann-Strafle 1, 85748 Garching, Germany 
^ Umversidad Complutense, Facultad de Fisicas, Ciudad Universitaria s/n, Madrid 28040, Spain 

Atomic quantum gases in the strong-correlation regime offer unique possibilities to explore 
a variety of many-body quantum phenomena. Reaching this regime has usually required 
both strong elastic and weak inelastic interactions, as the latter produce losses. We show 
that strong inelastic collisions can actually inhibit particle losses and drive a system into 
a strongly-correlated regime. Studying the dynamics of ultracold molecules in an optical 
lattice confined to one dimension, we show that the particle loss rate is reduced by a factor 
of 10. Adding a lattice along the one dimension increases the reduction to a factor of 2000. 
Our results open up the possibility to observe exotic quantum many-body phenomena with 
systems that suffer from strong inelastic collisions. 



Strong interactions are responsible for many interest- 
ing quantum phenomena in many-body systems: high- 
Tc superconductivity excitations with fractional 
statistic 0, topological quantum computation Q, and 
a plethora of exotic behaviors in magnetic systems [3]. 
One of the main physical mechanisms that gives rise to 
strong correlations for bosonic particles can be under- 
stood as follows. At low temperatures and for strong 
elastic repulsive interactions, particles tend to stay far 
away from each other in order to keep the energy low. 
That is, the wavefunction describing the particles tends 
to vanish when two of them coincide at the same position. 
In order to fulfill these constraints, this wavefunction has 
to be highly entangled at all times, which may give rise 
to counter-intuitive effects both in the equilibrium prop- 
erties as well as in the dynamics. In ID, for example, 
this occurs in the so-called Tonks-Girardeau gas (TGG) 
0) 1 where the set of allowed wavefunctions for bosonic 
particles coincide (up to some transformation) with those 
of free fermions. Despite being bosons, the excitation 
spectrum, the evolution of the density distribution, etc, 
correspond to those of fermionic particles. In 2D the 
same mechanism leads to the fractional quantum Hall ef- 
fect f?\ , where the ground state as well as the low energy 
excitations fulfill the above mentioned constraint, giving 
rise to the existence of anyons which behave neither like 
bosons nor fermions, but have fractional statistics [2]. 

Here we show that inelastic interactions can be used 
to reach the strong correlation regime with bosonic parti- 
cles: This may seem surprising because inelastic collisions 
are generally associated with particle losses. This behav- 
ior can be understood by using an analogy in classical 
optics, where light absorption is expressed by an imagi- 
nary part of the refractive index. If an electromagnetic 
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wave impinges perpendicularly on a surface between two 
media with complex refractive indices ni and n2, then 
a fraction \{ni — n2)/(ni -|- n2)P of the intensity will be 
refiected. In the limit \n2\ — * oo, the light is perfectly 
reflected off the surface, irrespective of whether n2 is real 
or complex. In our case, bosons interacting with large 
imaginary [1, [9| scattering length almost perfectly reflect 
off each other for an analogous reason, thereby giving rise 
to the same constraints in the particles' wavefunction as 
the ones corresponding to elastic collisions, and thus to 
the same physical phenomena. In our experiment, the 
correlations manifest themselves in a strong suppression 
of the rates at which particles are lost due to inelastic 
collisions. 

Our experiment uses molecules confined to ID by an 
optical lattice, both with and without a periodic poten- 
tial along the ID. We start with the transfer of a BEG of 
^^Rb atoms into a 3D optical lattice, in such a way that 
the central region of the resulting Mott insulator contains 
exactly two atoms at each lattice site. A Feshbach reso- 
nance at 1007.4 G [l3| is used to associate the atom pairs 
to molecules [ll| . Subsequently, the magnetic field is held 
at 1005.5 G. Atoms remaining after the association are re- 
moved with blast light. This procedure prepares a quan- 
tum state that contains one molecule at each site of a 3D 
optical lattice [H, [3 ■ The optical-lattice potential seen 
by a molecule is —V± cos^(A:a;)— VI cos^(fcy)— V|| cos^(A:2;), 
where A = 27r//c — 830.440 nm is the light wavelength. 
At the end of the state preparation V\\ ~ V± = 127 Er, 
where Er = h'^k'^/2m is the molecular recoil energy and 
m is the mass of one molecule. 

After state preparation, Vn is linearly ramped down to 
its final value. After this ramp, we have an array of tubes 
of ID gases. We choose a ramp duration of 0.5 ms. For 
much faster ramps, we observe a substantial broadening 
of the momentum distribution along the tubes. For much 
slower ramps, particle loss during the ramp becomes no- 
ticeable. After the ramp, the system is allowed to evolve 
for a variable hold time at the final value of V|| . During 
this hold time, molecules collide inelastically, leading to 
loss. After the hold time, all molecules are dissociated 
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FIG. 1: Time-resolved loss of molecules at Vji =0. The 
loss begins ai t — 0. The solid line shows a fit of Eq. [2] to 
the experimental data (•) with t < 1 ms. The best-fit value 
is x^^(O) = 4.3/ms, corresponding to Ksd = 2.2 x 10"^" 
cm^/s and, at t = 0, to g^^^ = 0.11. The dashed line shows 
the expectation for an uncorrelated system. The observed 
loss is much slower than the dashed line because of strong 
correlations. 



[ll| into atom pairs using the Feshbach resonance. The 
dissociation terminates the loss. Finally, the magnetic 
field and the lattice light are switched off simultaneously, 
and the number of atoms is determined from a time-of- 
flight absorption image. 

The experimental data represent an average over a 
large number of tubes of different lengths. This is not 
critical because the initial ID density n(0) = 2/A is iden- 
tical in all tubes. The lattice beams that create V± have 
a finite waist. This results in a harmonic confinement 
with angular frequency o^n = 27r x 71 Hz along the tubes. 
This is negligible as long as we evaluate the loss only 
for much faster timescales. In the decay of the molecule 
number as a function of hold time at =0 (Fig. [1]) 
the ramp down of begins at t — —0.5 ms and ends 
at t = 0. The data do not show noticeable loss during 
the ramp down. In order to avoid complications due to 
the harmonic confinement lo^^ along the tubes, we process 
only data for t < 1 ms. 

A quantitative understanding of the loss process is 
based on the ID particle density n, which evolves ac- 
cording to [T^ 



dn 
It 



(1) 



where g^"^") — (n^)/(n)^ is the pair correlation function 
that gives the reduction factor of the loss rate compared 
to an uncorrelated state where g'^^'' = 1. K is the ID 
loss rate coefficient, which can be related to the 3D scat- 
tering properties as follows: First, the scattering poten- 
tial can be modeled as a delta interaction with ID in- 
teraction strength g, yielding K — —21m(q)/h [13] ■ Sec- 
ond, extending the arguments of Ref. [l^ to the case 
of a complex-valued 3D scattering length a, one obtains 
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FIG. 2: Loss at V|| / 0. Solid hues show fits of Eq. |3]to the 
experimental data (filled symbols). Open symbols show re- 
sults of our numerical calculations. Black squares, red circles, 
and green triangles correspond to Vji/iJr = 1.8, 3.9, and 6.0 
respectively. 



g — 2?i^a/(ma^[l + aC (^) /v^ a±]), where a± is the size 
of the gas in the perpendicular direction and ( denotes 
the Riemann zeta function with C(i) « —1.46. The real 
and imaginary parts of a represent elastic and inelastic 
scattering, respectively. 

For the ground state of the TGG, one can derive an an- 
alytic expression for g^^K Introducing a dimensionless in- 
teraction strength 7 = mq/fi^n and using the same tech- 
niques as in Refs. [1^, [l3| one obtains g^'^^ = 47r^/3|7p 
in the limit I7I 3> 1. Inserting this expression in Eq. [T] 
we obtain 



dn 



(2) 



where x = 47r ?i K/5m \g\ . This equation gives a pre- 
diction for the particle losses for a low energy TGG. Simi- 
lar considerations apply to inelastic three-body collisions 
[il,[ii,i2fl]. Note that a TGG at finite temperature has 
recently been observed experimentally [111, |22, [1^ but for 
the known case of strong elastic interactions. 

A fit of Eq.[2]to the data in Fig. [T] with t < 1 ms yields 
Im(l/a) = 1/(24 nm). Combined with |Re(a)| < |Im(a)| 
[13] this yields 5(2) = g.ll ± 0.01 at t = 0. This shows 
that the loss rate is strongly reduced due to correlations. 
5'^^^ oc n^ decreases with time, so that the loss rate is 
reduced even further. 

From these results we can determine the two-body 
loss coefficient for a Bose-Einstein condensate K^d = 
-8nhlrR{a)/m = (2.2 ± 0.2) x IQ-^" cm^/s. This agrees 
fairly well with the measurements discussed further be- 
low and with a previous measurement [13] at 1005.8 G 
that yielded K^u — 1.5 x 10"^*^ cm^/s. For comparison 
we note, that if losses were not inhibited, i.e., g*^^^ — 1, 
then the loss should follow the dashed line, which is cal- 
culated with K^u from Ref. [13] ■ Clearly, the loss rate is 
reduced. 

This reduction is directly related to the fact that the 
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FIG. 4: Particle losses in a ID lattice. The initial level |1) 
contains exactly one particle at each lattice site. State |2) is 
obtained after one tunneling process. Population in state |2) 
decays incoherently into state |3) with rate P. The tunneling 
coupling between states |1) and |2) can be described by a 
Rabi frequency Q and a detuning A. In the limit P 2> 
the effective decay rate from |1) |3) is Peg ~ f2^/(P[l + 
(2A/P)^]) and the population of level |2) after a time t > 1/P 
is Pcft/P times that of level |1) [2^. The parameters are 
related to the Bose-Hubbard parameters by P = — 2Im(l7)/?i, 
A = Re{U)/h, and Q = VSJ/h. 



FIG. 3: Loss at different lattice depths Vji . A Fits as in Fig. 
[2]yield the experimental results (•) for n. A fit of Eq. |4]to the 
data yields the solid line. The best-fit value is Ksd = 1.7 x 
fO~^° cm^/s. The experimental data and the analytic model 
agree well with results of our numerical calculations (o). For 
comparison, the dashed line shows the naive estimate 2J/h 
which is nowhere even close to the data. B Pair correlation 
function g^^' = k/F calculated from the data in A. 



spatial wavefunctions of the particles do not overlap, and 
thus they become strongly correlated. Further support 
for this conclusion comes from the time dependence of 
our data. If the system were weakly correlated, then 
(/(^) w 1 would be time independent. Eq.[T] would then 
predict that the number of particles N follows dN/dt cx 
instead of dN/dt cx N'^ if we are close to the TGG 
ground state • We fit dN/dt oc N^ with an arbitrary 
power p to the data with t < 1 ms. This yields p = 
4.3 ± 0.6 in good agreement with p = 4. 

We now turn to a situation with a lattice potential 
along the ID tubes with V\\ ^ VI. The motion perpen- 
dicular to the tubes remains frozen out as before, but 
the motion along the tubes is now described as hop- 
ping between discrete lattice sites. This is a natural 
way of amplifying the effects due to interactions and thus 
reaching more deeply into the strong-correlation regime 
(2lll25|. Apart from that, many paradigmatic models in 
solid state and other fields of physics assume a lattice 
structure. We performed measurements (Fig. ^ similar 
to that in Fig. [1] for various values of the lattice depth 
V||. Here, the temporal change of g'^'^'^ is negligible 
Spatial integration of Eq. [1] yields 



where we abbreviated k = Kn{0)g'^^\ We fit this to the 
data with N{t) > N{0)/2 in order to neglect the har- 
monic confinement along the tube. The best-fit values 
for the loss rate k obtained for different lattice depths V|| 
are shown in Fig. [3lA_. 

A prediction for n can be developed starting from the 
Lieb-Liniger Hamiltonian [ij, 13 • In the presence of 
a periodic potential in the tight-binding limit, one ob- 
tains a Bose-Hubbard model [26] with tunneling ampli- 
tude J and on-site interaction Ke{U). J and U can be 
expressed in terms of m, a, and a Wannier function, 
r = —2lm{U)/fi is the rate at which two particles at 
the same site are lost. In this model losses occur when 
two particles occupy neighboring sites and one of them 
hops to the other's site. For strong inelastic interactions 
{J/hT ^ 1) one obtains an effective loss rate Fcff for 
two neighboring particles (Fig. [4|) . An extension of this 
double-well model to many sites yields k = 4reff so 
that 
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Im(a) 



(4) 
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iV(0) 



N'^{t) 



(3) 



This coefficient shows again that for strong inelastic in- 
teractions (r — > oo) particle losses are inhibited. For 
large, finite values of F one can, as before, observe such 
an inhibition by looking at the particle losses. We fit 
Eq. Hto the data in Fig. d^. With |Re(a)| < |Im(a)| 
there is only one free fit parameter. The best-fit value is 
K:iD = (1-7 ± 0.3) X 10^1° cm^/s, which is close to the 
result of Fig. [TJ As V\\/Er increases from 1.7 to 10, F 
increases from 45/ms to 82/ms. 

In order to test the quality of this analytic model, 
we performed extensive numerical calculations by solv- 
ing the master equation [T^ using matrix product den- 
sity operators [27| . They reveal that the system maps 
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to a very good approximation to a fermionized gas and 
that it loses its memory about the initial state in a time 
^ 1/r. During this very short time only little loss occurs, 
and after this short transient the loss is well described by 
a time independent 5^^^ with k bom Eq. |4l We find good 
agreement between experimental data, analytic model, 
and numerical results in Figs. [2] and [H 

For |Re(a)| < |Im(a)|, Eq.[3]tells us that the larger the 
loss coefficient F, the smaller the actual loss rate k. This 
means that fast on-site loss tends to preserve the initial 
state and thus suppresses tunneling in the many-body 
system. This can be interpreted as a manifestation of 
the continuous quantum Zeno effect [281] : fast dissipation 
freezes the system in its initial state. Without this Zeno 
effect, one might naively estimate that tunneling would 
occur at a rate ~ 2J/h. If each such tunneling event 
would lead to immediate loss, then 2J/'h should set the 
timescale for the loss, but that estimate is too nai've, as 
the dashed line in Fig. [3K shows. 

The value of K^d extracted from Fig. is used to 
calculate F and thus g^^^ = k/T for each exper- 
imental data point. The results (Fig. [3j3) agree well 
with the theoretical expectation (solid line) based on 



the same value of K^d. The smallest measured value 
of g^'^'> = (4.6 ±0.7) X 10^'' represents an improvement of 
more than two orders of magnitude over previous exper- 
iments 23]. Note that a noninteracting gas has g^^^ = 1 
at any lattice depth. The observed suppression of g^^^ is 
caused by the inelastic interactions, not by the lattice it- 
self. Moreover, the strength of the correlations in our ex- 
periment is determined by the interparticle interactions. 
This differs fundamentally from experiments in very deep 
lattices where tunneling is negligible on the timescale of 
the experiment, interactions are irrelevant, and bosonic 
symmetrization of the wave function is possible but has 
no detectable consequences. 

The mechanism introduced here could also give rise 
to other strongly-correlated states, such as a Laugh- 
lin state ,7,] or one with anyonic excitation 2]. The 
present paper opens up the possibility of observing ex- 
otic quantum many-body phenomena in systems that 
suffer from strong inelastic collisions. Furthermore, the 
rate coefficients for those collisions may be artificially in- 
creased using photoassociation or Feshbach resonances, 
thus further reducing the actual loss rate in the strongly- 
correlated regime. 
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MODELING DETAILS 

Real Part of the Scattering Length 

First, we analyze the contributions of Rc(a) and Im(a) 
to the value Im(l/a) = 1/(24 nm) extracted from Fig. [T] 
The identity 



Re(a)2 = -Im(a)2 - 



Ini(a) 
Im(l/a) 



(SI) 



shows that Re(a)^ is a parabola as a function of Im(a) 
at fixed Im(l/a). This yields an upper bound |Re(a)| < 
l/(2Ini(l/a)) = 12 nm. This information can be com- 
bined with our previous studies of excitation spectra 
[l^ . [l^. With the upper bound for |Re(a)|, the exci- 
tation spectra should yield a resonance at a frequency 
below 1.4 kHz for V\\ = V_i_ = 15 i?,.. We performed 
measurements that cover a broad range of experimental 
parameters but could not detect such a resonance. This 
can only be explained if |Re(a)| <^ |Im(a)|, because in 
this case the resonance becomes very broad and shal- 
low, so that it might be undetectable. We conclude that 
|Re(a)| < |Im(a)|. As a result Im(l/a) w -1/Im(a). 

The ID measurements (no matter if V|| =0 or 7^ 0) 
yield an experimental value for Im(l/a), while the 3D 
measurement in Ref. [13] yields a value for Im(a). In 
principle, these two values could be inserted into Eq. EH 
to extract a value (not just an upper bound) for |Rc(a)|. 
But in practice the systematic uncertainty in Ref. [24] is 
so large, that we cannot constrain |Re(a)| any further. 



Dissipative Lieb— Liniger Model 



Now we turn to modeling the loss in the ID tubes at 
V|| = 0. A ID spinless bosonic system in a box with pe- 
riodic boundary conditions (similar conclusions apply to 
many other situations) is described by the Lieb-Liniger 
Hamiltonian [la| 



2m ^ dxf 



(S2) 



where Xi is the position of the ith boson. Generaliz- 
ing the techniques developed in Ref. [lB| to the case of 
complex-valued g, one can show that the imaginary parts 



of all eigenvalues of H vanish as \g\ 00 and that the 
corresponding eigenfunctions no longer overlap. The for- 
mer indicates that in this limit there will be no losses, 
while the latter gives rise to strong-correlation phenom- 
ena. For small values of the relative position of two par- 



ticles, 



the eigenfunctions have to scale 
2 



as -0 oc 2h^/mg + \xij \ + 0{x'jj) in order to satisfy the 
Schrodinger equation. In the above limit and for Xij — 0, 
■0 must vanish regardless of the argument ol g. In any 
experimental situation, \g\ will be finite. In that case, 
there will be finite losses but they will remain small if |g| 
is large. The master equation discussed further below can 
be used to derive the equation describing the evolution 
of the ID particle density n 



d{n) 
dt 



-Kg 



(S3) 



where K = —2lm{g)/Ti is a ID loss rate coefficient and 
_ ^^2y ^^-^2 |-]-^g pg^jj. correlation function. Time 

integration of Eq. [2] followed by spatial integration yields 

N{t) = N{0){1 + 3n3(0)xt)"^/^. 



Dissipative Lattice Model 

We now turn to a situation with a lattice along the 
ID tubes and with periodic boundary conditions. The 
time evolution of the particle number Uk at site k can 
be derived from the master equation discussed further 
below. The result closely resembles Eq. [S3l 



djrik) 
dt 



(S4) 



where in the lattice g^^^ — {nk{nf; — l))/(nfe)^. In or- 
der to obtain an analytic approximation for 5'^' we con- 
sider the limit J <ti hV, where the probability to have 
more than two particles at one site is negligible. Hence 
gC^) = 2p2/p\ where pi is the probability of having i 
particles at site k. We then use the three-level model 
of Fig. m We can write p2 = gnFeff/r, where qu is 
the probability that both site k and one of its neighbors 
is occupied, which we can estimate as gn ~ 2p\. This 
approximation is motivated by the observation that for 
J <C TiF the three-level model of Fig. g] yields hVcS < J- 
Hence, the particles have time to redistribute across the 
lattice between successive loss events. Our numerical re- 
sults confirm that this is a good approximation in the 
parameter regime of our experiment. Thus we get 



9 



(2) 



k/F 



(S5) 



6 



with K = 4rcff. The approximation qn ~ 2pl makes 
f/f^^^ independent of density and therefore independent of 
time, in contrast to the case V|| = 0. g^"^^ can be made ex- 
tremely small by reducing the tunneling amplitude, and 
thus one can approach a fermionized gas. Furthermore, 
assuming that each lattice site is initially occupied by one 
particle, spatial integration of Eg. IS4I yields 



dN 

IF 



iV(0) 



N'^{t) 



(S6) 



Time integration yields N(t) = ^(0)/(l + Kt). 



Master Equation 

In general, particle losses are described in terms of a 
master equation for the density operator p describing the 
particles 



n'^ = -iHp + ipH^ + J[p) 
at 



(S7) 



For the case without the lattice along the ID tubes, H is 
the second-quantized version of Eq. IS2I and 



J{p) = -Im(g) / dxi/^{x)p'^''^{x) , 



(S8) 



where ^{x) is the bosonic field operator. If we assume 
|Im((7)| — > oo, then the wavefunctions do not overlap and 
thus i/(jo) = 0. Thus, the master equation is reduced 
to a simple Schrodinger equation (for an initially pure 
state), hd\'i>)/dt — —iH\'i'), where we can make use of 
the fermionization transformation of Girardeau. On the 
other hand, for finite losses, we can make use of the mas- 
ter equation to derive the evolution of the total particle 
number 



N ^ dx{'ii''{x)i'{x)) 



(S9) 



obtaining Eq. IS3I for a homogeneous system. 

In the presence of the lattice along the ID tubes, we 
obtain a Bose-Hubbard model 



H = -jY^{alak+i+H.c.) + -Y^afal (SIO) 

k k 

J{p) = -Im([/) ^a^paf (Sll) 



Finite Energy in the Experiment 

Finally, we address the question, how close our exper- 
iment is to the ground state of the TGG. Eq. [2] assumes 
that we have a TGG at low energy, close to the ground 
state. Thus, the reasonable agreement of our experi- 
mental data with that equation gives a strong indication 
that indeed we are relatively close to the ground state. 
As the loss proceeds, the particle number changes and 
the ground state changes correspondingly. A simple es- 
timate for the temporal evolution of the energy can be 
obtained when assuming that the kinetic energy per par- 
ticle E^™ /N would remain constant when two particles 
are lost. This is to be compared to the new energy of 
the ground state which now has two less particles. In 
a harmonic oscillator, the kinetic energy of the fermion- 
ized ground state scales as E^™ (x N'^. As the loss pro- 
ceeds, the system thus automatically evolves away from 
the ground state, because E^'''{t)/Ef''{t) = N{0)/N{t). 
Actually, this gives us an upper bound to the energy. The 
reason is that in the limit studied here, for the (pseudo) 
eigenstates of Eq. [S2] the imaginary part of the corre- 
sponding eigenvalues is proportional to the real (kinetic) 
energy. Thus, the states with higher kinetic energies are 
preferentially depleted. In any case, ignoring this last 
fact we estimate for the data with t < 1 ms in Fig. [1] that 
^kin^^kin jjjcreases by a factor of ~ 2 during the loss. 
This is comparable to the energy of the TGG observed in 
Ref. [2l[, where temperatures of the order of the Fermi 
energy were achieved. 



and the same as before can be done. 

Our numerical calculations are based on this master 
equation. As J ^ hT, they neglect the possibility that 
more than two particles could occupy the same site. The 
lattice depth is modeled as time-independent and the 
evolution of the particle number starting from a quantum 
state with one particle at each lattice site is modeled 
without any assumptions about the time dependence of 



